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Applying the 1 + 3 formalism we write down the full set of equations governing the structure and 
the evolution of self-gravitating cylindrically symmetric dissipative fluids with anisotropic stresses, 
in terms of scalar quantities obtained from the orthogonal splitting of the Riemann tensor (structure 
scalars), in the context of general relativity. These scalars which have been shown previously (in 
the spherically symmetric case) to be related to fundamental properties of the fluid distribution, 
such as: energy density, energy density inhomogeneity, local anisotropy of pressure, dissipative flux, 
active gravitational mass etc, are shown here to play also a very important role in the dynamics 
of cylindrically symmetric fluids. It is also shown that in the static case, all possible solutions to 
Einstein equations may be expressed explicitly through three of these scalars. 
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I. INTRODUCTION 



In a recent series of papers [lH4( a set of scalar func- 
tions obtained from the orthogonal splitting of the Rie- 
mann tensor and referred to as structure scalars, were 
introduced in the discussion about the structure and 
evolution of spherically symmetric fluid distributions. 
Such scalars (five in the spherically symmetric case) were 
shown to be endowed with distinct physical meaning. 

In particular they control inhomogeneities in the en- 
ergy density [l[ , and the evolution of the expansion scalar 
and the shear tensor [IH1. Also in the static case all pos- 
sible anisotropic solutions are determined by two struc- 
ture scalars 

Furthermore, the role of electric charge and cosmolog- 
ical constant in structure scalars has also been recently 
investigated 

It is our purpose in this work to carry on a study 
on cylindrically symmetric fluid distributions based on 
structure scalars. The motivation to undertake such an 
endeavour is provided on the one hand by the conspic- 
uous physical meaning of such scalar quantities, and on 
the other hand by the instrinscc interest of cylindrically 
symmetric systems in general relativity (see and 
references therein). 

For doing so we shall apply the 1 + 3 formalism devel- 
opped in [2(>t24j . However we shall not follow a frame 
formalism but a coordinate basis approach in which the 
orthonormal frame is only used to identify frame compo- 



nents of proper vectors as scalars that can have a covari- 
ant interpretation. 

We shall first define the structure scalars corresponding 
to a general cylindrically symmetric fluid distribution, 
then we shall deploy the full set of equations governing 
the structure and evolution of such a system and express 
them in terms of the above mentioned scalars. A sys- 
tematic though non exhaustive study of these equations 
is carried out, including the coupling of the generalized 
"Euler" equation with a transport equation. 

Besides the structure scalars, we shall also introduce 
a set of scalars describing the shear tensor and the mag- 
netic and electric parts of Wcyl tensor. 

A subset of our system of equations will be shown to 
exhibit the role of structure scalars on the evolution of 
the expansion scalar and the shear tensor. Another sub- 
set will relate the shearfree condition, the dissipative flux 
and the magnetic part of the Weyl tensor, while the last 
four equations of our system determine the inhomogene- 
ity factor and its evolution. 

The static case shall be considered in detail. In general 
it will be shown that any static solution is defined by a 
triplet of structure scalars. 

Finally all results will be discussed in some detail in 
the last section, and a list of issues deserving further 
attention will be presented. 



II. COLLAPSING ANISOTROPIC DISSIPATIVE 
FLUID CYLINDERS: BASIC DEFINITIONS AND 
NOTATION 
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We shall start by introducing basic definitions and no- 
tation to be used throughout the text. We shall closely 
follow the notation of [Tj] except for the fact that we are 
now considering a dissipative fluid. We assume the gen- 
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eral time dependent diagonal non rotating cylindrically 
symmetric metric 



ds 2 



-A 2 (dt 2 - dr 2 ) + B 2 dz 2 + C 2 d(j> 2 , (1) 



where A, B and C arc functions of t and r. To 
represent cylindrical symmetry, we impose the following 
ranges on the coordinates 

-co < t < oo, < r, -co < z < co, < cj) <2n, (2) 

where we assume C = at r = which is a non-singular 
axis. We number the coordinates x° = t, x 1 = r, x 2 = z 
and x 3 = (p. 

Next, let us consider a collapsing cylinder filled with 
anisotropic and dissipative fluid. Our study concerns 
cither bounded or unbounded (cosmological) configura- 
tions, in the former case we should further assume that 
the fluid is bounded by a timelike cylindrical surface E. 

Thus, the energy momentum tensor is given by 

T a p = Qi + P)V a Vp + Pg af 3 + q a V fj + qpV a + U aP , (3) 
where 



Now, in order to avoid misunderstandings, the follow- 
ing remark is in order: It is well known that the choice of 
the vector V a is not unique, we could for example choose 
a "tilted" congruence and therefore the splitting of the 
energy-momentum tensor would be different of course. 
In this sense our study is related to the congruence of 
observers which are at rest at each point with respect 
to the corresponding fluid clement, i.e. a congruence for 
which the four velocity of the fluid is V a given by ([9]). 

With the notation above we can write nonvanishing 
components of the Einstein equations (see Appendix I). 



A. Kinematical Variables 

Since we are considering nonrotating fluid distribu- 
tions, there are only three kinematical variables: the ex- 
pansion O, the four acceleration a a and the shear a a a 
which as usual are defined by: 



6 = V c 



(10) 



u a0 = U s {S a S p - h aP ) + U k (K a K p - ^h a0 ), (4) 



a a = V a ., V . 



(11) 



h a p = g a p + V a V/3, 



q a — qL ai 



(5) 



(6) 



and scalars II S , Ilfc and q, are functions of t and r. 

Alternatively, we may write the energy momentum ten- 
sor in the form 

T a p = {n + P r )V a Vp + P T g a p + q{L a V, 3 + LpV a ) 

+ (Pz - Pr)S a S (j + (P - P r )K a Kp, (7) 

where 



( Pz - Pr 



IT 



(P* 



-P r ) = U k , 

Pz + Pr 



(8) 



The unitary vectors V a ,L a ,S a , K a form a canonical 
orthonormal tetrad. V a is a hypersurface orthogonal 
4— velocity vector, S a ,K a are tangent to the orbits of 
the 2— dimensional group that defines cylindrical sym- 
metry and L a is orthogonal to these orbits and to V a . 
With the above definitions it is clear that ^ is the energy 
density (the eigenvalue of T a p for eigenvector V a ), q a is 
the radial heat flux, whereas P is the isotropic pressure. 

We choose the fluid to be comoving in this coordinate 
system, hence 



V a = -AS° a , L a = A6l S a 



B6Z, K a 



CSl (9) 



c Q/ 3 = V( Q;(3) + a( a Vs) - -Oh al3 . (12) 



Using ((T|), (fTU|) and (fT2"]) we obtain for the expansion, 

(13) 



Q 1 A B C 

q= a\a + b + c 



and for the non zero components of the shear, 



°33 
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(14) 
(15) 
(16) 



where overdots and primes stand for partial differentia- 
tion with respect to t and r respectively. 
We can also express the shear tensor as 



S a S(S — —h a p J + (7k (K a Kp — —h a p 



with 




(17) 
(18) 
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and 




CTsCTfc 



(19) 



(20) 



Observe that unlike the spherically symmetric case, the 
shear depends now on two nonvanishing independent 
scalars. Finally, for the four acceleration a a we obtain 



A' 

a„ = aL„ with a = —rr. 

A 2 



(21) 



Before ending this section it should be recalled that in 
the case of bounded configurations, junction conditions 
(Darmois) should be satisfied on the boundary surface in 
order to exclude the presence of thin shells. Such con- 
ditions shall not be deployed here since we shall not use 
them but the reader may found a detailed discussion on 
this point in [ljj, [l7|. 



III. THE ORTHOGONAL SPLITTING OF 
RIEMANN TENSOR AND STRUCTURE 
SCALARS 

In order to introduce the structure scalars correspond- 
ing to our problem, let us first recall that the Riemann 
tensor may be expressed trough the Weyl tensor C p a/3Al , 
the Ricci tensor R a p, and the scalar curvature R. 

The electric (E a p) and magnetic (H a p) parts of the 
Weyl tensor C a p 7 Si are defined as usual by 



E a p — CavpsVV 5 , 



(22) 



where rj auep is the Levi-Civita tensor and e a pp = 
Calculating (|22|) for the metric ([l} we obtain 

E a p = E s (S a Sp — -hep) + E k (K a Kp — -h a p), (23) 



with 



Es = A r B 2 ~ C ° 202 ~ A^ C ° 1QU 

Ek = JpQ2 C ° 303 ~ ^4^0101, (24) 



and 



with 



H a p = H{S a K p + SpK a ), 



H 



a 



0313 



A 2 C 2 



(25) 



(26) 



where the explicit expressions for the Weyl tensor com- 
ponents (computed for the orthonormal "canonical tetr- 
dad") may be found in the Appendix II. 

It is worth observing that now, unlike the spherically 
symetric case, the magnetic part of the Weyl tensor does 
not vanish (in general) and the electric part depends on 
two independent scalar functions. 

Now, the orthogonal splitting of the Riemann tensor is 
carried out by means of three tensors Y a p , X a p and Z a p 
defined as (see [3, [25|, [2(| for details) 



and 



X a p = -jV a ^ PR * P ps vUvt \ 



Z a p — TfatpRsp^V 1 ', 



(27) 



(28) 



(29) 



where R* a p„ s = lr] epvS R a p £f '. 

Using the decomposition of the Riemann tensor in 
terms of the matter variables and the electric and mag- 
netic parts of Weyl tensor (seeQ, H3) and ((23]) ([27]>. we 
obtain: 

Yap = -^Yrh-ap + Y s (S a Sp — -h a p) 

+ Y k (K a Kp-^h a p), (30) 



with 



Yt = -(n + P z +P<f> + Pr), (31) 
Y S = E S -^{P Z -P r ), (32) 
Y k = E k -~{Pt-P r )- (33) 
In a similar way the tensor X a p can be written as: 

X a p = -^Xxhap + X s (S a Sp — -h a p) 

+ X k (K a Kp-h a p), (34) 



with 



Xt — Kfi, 

X s = -E s -^(P Z -P r ), 

X k = -E k -~(Pt-Pr). 

Finally, from ([22]). {25) and $23$ we obtain 
Z a p = H a p + -^Kq p e a p P - 



(35) 
(36) 

(37) 



(38) 
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From (|38|) two scalar functions may be defined as fol- 
lows: 



Z H =2H = (S a K p + S p K a )Z, 



a/3, 



(39) 



one from the dissipativc flux and the other from a term 
which represents gravitational radiation [2^, [2i| , this last 
term being proportional to the magnetic part of the Weyl 
tensor. 



nq 



(S^K* - S a K f) )Z a[: 



(40) 



Thus the full set of structure scalars are now defined 
by the eight quantities: Y T , Y s , Y k , X T , X s , X k , Z H , Z q . 
The corresponding expressions of these scalars in terms 
of the metric functions are given in the Appendix III. 

Before ending this section it would be useful to in- 
troduce a relevant quantity defined in terms of tensors 
Y a p, X a j3, Z a p. This is the super-Poynting vector defined 
by 



(y7 Z /M _ X ]Z 5 P) , 



(41) 



in our case the above expression becomes 

.2 



P a = 2H(E k - E s )S p K' 1 e a p 1 + y (ji + P r )q a . (42) 



IV. BASIC EQUATIONS 

In this section we shall deploy the relevant equa- 
tions for describing a dissipativc self-gravitating locally 
anisotropic, cylindrically symetric fluid. In spite of the 
fact that not all these equations are independent (for ex- 
ample the field equations and the conservation equations 
(Bianchi identities)) we shall present them all, since de- 
pending on the problem under consideration, it may be 
more advantageous using one subset instead of the other. 
As mentioned in the Introduction these equation are ob- 
tained applying the 1+3 formalism to cylindrical symme- 
try |20l - [24| (for the specific case of spherical symmetry see 



A. Conservation Laws 



This four-vector describes the flux of superenergy and, 
as it is apparent from (|42[) . it embodies two contributions: 

I 



The conservation law T a Jr = leads to the following 
couple of equations 



M +-(P r +P z +P ) 



Wq + q 



2a 



1 ( B' C 



C 



A \ B 



i(P z - P r )(2a 3 - a k ) + \{P<t>- Pr){2a k - <r.) = 0, 



(43) 



VP r + D t q - 



A 



(Pz - Pr 



B' 



Pr) 



c 



+ (jj. + P r )a - -(a s + cr fe - 49)g = 0, 



(44) 



r 



where D t f = f, a V a and V/ = /, Q L C 



equation) and the shear tensor, as well as some constraint 
equations are obtained: 



B. Ricci Identities 

From the Ricci identities for the vector V a the following 
evolution equations for the expansion (the Raychaudhuri 

I 



D t & -Va-a 2 -a±f^ + ^pj+ ±Q 2 + ^(a 2 s - a s a k + a 2 k ) = -Y T , 



(45) 



a j [2^ -^L\- a 2 -Va- D t (2a s - a k ) + ^(a 2 - 2a 2 + 2a s a k ) - ^Q(2a s - a k ) 



2Y S -Y k , (46) 



r 

J 



1 ( B 1 C'\ 1 
2Va + 2a 2 - a- f — + — J + D t (a s + a k ) + -(a 2 + a\ - 4a s a k ) 



~Q(a s +a k ) = -(Y a + Y k ), (47) 



1 B' C 1 

-V (26 + * a + * k ) +cr s — + o k — = Z q , (48) 



C. Differential equations for the Weyl tensor 



C B' 



The Bianchi identities, using the Einstein equations, lead 



Z H = V(cr fe - a a ) + ( Vk-^j - G »~^j > ( 49 ) to tnc following equations: 
I 



V(y s +Y k -X s - X k ) - 3(Y S - X s )^ - 3(Y k - X k )^ - 6H(a s - a k ) 
= «V(2jtx + P r + P z + P ) + 3k(^ + P r )a + 2nq(Q - a s - <j k ) + 3kD ( q 



(50) 



r 

j 



V(2F S -Y k - 2X S + X k ) + 3(Y S ~X a )^ + 3a(Y s -Y k - X s + X k ) + 6H{G - a k ) 
+6D t H = -kV(/x -P r -P z + 2P <t> ) - 3k(P - Pr)jg + w{® - a, + 2cr fe ), 



(51) 



r 

j 



a s (-Y s + 2Y k + X s - 2X k ) + a k (2Y s -Y k - 2X S + X k ) + 6(Y S + Y k - X s - X k ) 

B' C 
AB ~ AC 

k(h + P r )(@ - (J s - o-fc) + kA(/i + 2P r -P z - P<p) + 3n\7q + 6nqa 



+D t (Y s + Y k -X a -X k )-6H 



(52) 
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a s (2Y s -Y k - 2X S + X k ) - a k (Y s + Y k - X s - X k ) - 0(2Y S -Y k - 2X S + X k ) 

-D t (2Y s -Y k - 2X S + X k ) - 6H - 6VH - \2Ha 



K(p + P Z )(Q + 2a s - a k ) + K D t {n - P r + 2P Z - P ) + 3nq 



B' 
AB' 



(53) 



r 



We shall next proceed to analyze different problems 
using different subsets of the equations above. 



V. DYNAMICAL EQUATION, TRANSPORT 
EQUATION AND THERMOINERTIAL EFECT 

We shall now elaborate on (|44p (the generalized "Eu- 
ler" equation) as follows. 

By analogy with the spherically symmetric case let us 
first define the "velocity" 



U 



c 

A 



D t C, 



(54) 



then using f|126[) we get 

C 

D t U = a-r- nP r C 

C (B AB BC B'C 
~l¥\B~~^B + B~C~ BC 



AB' 
AB 



,(55) 



that can be also written down as (see Appendix I) 



C 

D t U = a—- kP t C 



R 



2323 
C 2 



(56) 



Solving the above equation for the a term and feeding 
this back into (l44l) we obtain 



A 



(ji + P r )D t U = -{n + Pr) 



pr C ( R.Q2V. 



R 



A 2 



2323 



X7P r + (P z -Pr)^ + (P* - Pr)^ 



A 



1 



-D t q+-(a s + a k -AQ)q 



(57) 



The equation above has the "Newtonian" form 



Force = Mass density x Acceleration. 

Indeed, the term on the left is the inertial mass (den- 
sity) multiplied by the proper time derivative of the "ve- 
locity" U. On the right hand we have three different 
terms: the first one is just the "gravitational force" term 
(see below), the second one includes hydrodynamic force 
terms (presure gradient plus the anisotropic contribu- 
tions), finally the last term represent the contribution 
from dissipative processes. 

Now, before proceeding further, let us identify the first 
term on the right of ([57]) as the "gravitational force" . For 
doing that, let us consider the static vacuum case (Levi- 
Civita). From Einstein equations f|124p — (|128[) it follows 
that the metric ([TJ in this case has components: 



p _ a 



C 



r (/3+l) 



a 



(58) 



(59) 



A 



r /3(/3+l) 



(60) 



with a, P constant. 
Then, the term 



C (R 



t-0202 -^2323 



becomes 



B 2 \ A 2 C 2 



T 2 ^ + l)2 - 



(61) 



(62) 



Now in the weak field limit the gravitational potential 
of an infinite line with mass per unit of length a is 



$ = 2a In r + constant. 
Therefore in that limit we have 
P 2 + P = 2a, 



(63) 



(64) 
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implying 



-1 + Vl + 8cr 



since in the weak field limit a << 1 then 
/3 k, 2a. 

Thus (|62p becomes in the weak field limit 
C 2a 



(65) 



(66) 



(67) 



but in that limit C ~ r, A sa 1 and the term above is 
just the gravitational force exerted by the infinite line. 



Thermoinertial effect 



Let us now get back to our equation (|57j) . In order to 
obtain an expression for the D t q term we have to resort to 
a transport equation. We shall need a transport equation 
derived from a causal dissipative theory ( e.g. the Miiller- 
Israel-Stewart second order phenomenological theory for 
dissipative fluids (3ll - [33 |V 

Indeed, as it is already well known the Maxwell- Fourier 
law for radiation flux leads to a parabolic equation (dif- 
fusion equation) which predicts propagation of pertur- 
bations with infinite speed (see [35|-[38( and references 
therein) . This simple fact is at the origin of the patholo- 
gies [3^1 found in the approaches of Eckart [40[ and Lan- 
dau |4l| for relativistic dissipative processes. To over- 
come such difficulties, various relativistic theories with 
non- vanishing relaxation times have been proposed in the 
past [3l|-BEllIl. Although the final word on this is- 
sue has not yet been said, the important point is that all 
these theories provide a heat transport equation which is 
not of Maxwell- Fourier type but of Cattaneo type 44j, 
leading thereby to a hyperbolic equation for the prop- 
agation of thermal perturbations (see [45| for a recent 
discussion on this issue). 



J 



Then feeding back (|69|) into (|57|) . we obtain 



a_ 

A 



VP r + {P z - P r ) 



B' 
AB 



ili + Pr) 



(P^ - Pr) 



KT 



A key quantity in these theories is the relaxation time 
t of the corresponding dissipative process. This positive- 
definite quantity has a distinct physical meaning, namely 
the time taken by the system to return spontaneously 
to the steady state (whether of thermodynamic equilib- 
rium or not) after it has been suddenly removed from 
it. Therefore, when studying transient regimes, i.e., the 
evolution from a steady-state situation to a new one, r 
cannot be neglected. 

Sometimes in the past it has been argued that dis- 
sipative processes with relaxation times comparable to 
the characteristic time of the system are out of the hy- 
drodynamic regime. However, the concept of hydrody- 
namic regime involves the ratio between the mean free 
path of fluid particles and the characteristic length of 
the system. Therefore that argument can be valid only 
if the particles making up the fluid are the same ones 
that transport the heat. However, this is never the case. 
Specifically, for a neutron star, r is of the order of the 
scattering time between electrons (which carry the heat) 
but this fact is not an obstacle (no matter how large the 
mean free path of these electrons may be) to consider 
the neutron star as formed by a Fermi fluid of degener- 
ate neutrons. The same is true for the second sound in 
superfluid Helium and solids, and for almost any ordinary 
fluid. In brief, the hydrodynamic regime refers to fluid 
particles that not necessarily (and as a matter of fact, 
almost never) transport the heat. Therefore large relax- 
ation times (large mean free paths of particles involved 
in heat transport) does not imply a departure from the 
hydrodynamic regime (this fact has been stressed before 
[4q |. but it is usually overlooked). 

Thus, the transport equation reads 



-Kh af} (Tp + Ta) 



-KT 



\KT 2 J 



(68) 



where K and T denote thermal conductivity and tem- 
perature respectively, and whose only component in our 
case is 



D t rq 



-K{VT + Ta) - l -KT 2 qD t 



-rqQ. (69) 



Then feeding back |69|) into (|57]l. we obtain 



t(h + P r ) 



D t U 



i^ + Pr) 



KT 



a 

AC 



A 



1 KT 2 ( t 
2— qDt \KT* 



T{H + P r ) 

KVT 1 



-P r C 3 



1 

C2 



: qO 



(70) 



where the function m is defined by Equation (|70[) indicates how inertial thermal effects re- 

fiP r C 3 C 3 ( i?0202 -R2323 \ r -~s 

m= —'m {-&—&-)' (71) 
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duce the effective inertial mass. This effect and its con- 
sequences, first reported in 1471] . has been extensively dis- 
cussed in the past (see [^M55l ] and references therein). 



VI. EVOLUTION OF THE EXPANSION 
SCALAR AND THE SHEAR: THE GEODESIC 
CASE 

Let us now turn to equations (|4"5"]) . ([41)]) and ([4T]) . The 
former is just the evolution equation for the expansion 
scalar (Raychaudhuri equation) for the cylindrically sym- 
metric case whereas the latter describes the evolution of 
the shear. For simplicity wc shall restrict ourselves to the 
geodesic case (a = 0). 

First of all observe that the evolution of the expan- 
sion scalar is fully controlled by the scalar Yt (as in the 
spherically symmetric case). For the shear however we 
have two equations (for the two independent components 
of the shear tensor). Introducing the variables: 



oi 



07c 



an = 2a s 



Ok, 



Y T = Y a +Y k ; Y n =2Y a -Y k , 
46]) and (147]) become (in the geodesic case) 



(72) 



(73) 



Dtoj = -Yi - i(2^j -oj- 2ai<r n ) - \&<ri (74) 



A 07/ 



-Yi 



ii 



laidn - (t'jj) - -Qou. 



(75) 



The two equations above describe the evolution of the 
shear tensor, which is fully controlled by the scalars Yj 
and Yjj (Y s ,Yk). If we assume that the fluid is initially 
(t = 0) shearfree and scalars Yj and Y// vanish for a time- 
like interval {t\ > t > 0), then it follows at once from ([74]) 
and ([75]) that the fluid remains shearfree in that inter- 
val. However, even small deviations from the vanishing 
condition of the above mentioned scalars would produce 
deviations from the shearfree condition. In other words, 
all the information about the stability of the shearfree 
condition is encoded in Yj and Yu (for a discussion on 
this problem in the spherically symmetric case see Q). 



VII. THE LINK BETWEEN SHEARFREE 
CONDITION, DISSIPATIVE FLUX AND THE 
MAGNETIC PART OF THE WEYL TENSOR 



Wc shall now try to extract some of the information 
contained in (gSJ) and (gSJ). 

First of all observe that ([4"5]) implies that under the 
shearfree condition the inhomogeneity of the expansion 
scalar V8 is controlled by Z q . Also from (|4l?|) it follows at 



once that the shearfree condition implies that the mag- 
netic part of the Weyl tensor vanishes, this last result 
was known for perfect fluids [5(| . 

Now, the remaining relevant question is: what can we 
infer about the shear from the vanishing of H and ql We 
were unable to elucidate this question in the general case, 
therefore in what follows, we shall restrict to the geodesic 
case (a = 0). Under this latter condition 



A! 

a = — = => A = A(t). 
A 2 v ' 



(76) 



Then assuming q a = H a p = 0, and taking into account 
76]) . we obtain from p25) and ([26]) 



~B 



c 



A 

A (~B 



H 



1 

2^2 



D 



C 



A 
A 



C 

&_ 

B 



9L 
c 



(77) 



0. (78) 



Combining and integrating the equations ([77]) and ([78]) 
we obtain 



B(t,r) = A(t)b(r) + a(t), 
C(t,r) = A(t)c(r) + p(t), 



(79) 
(80) 



where b, c, a, (3 are arbitrary functions of their argument. 

Next, from the regularity condition C(t, 0) = and 
redefining the function c(r) we may write 



C(t,r) =A(t)c(r). 



(81) 



From the above equation and ([T9"]) it follows that = 
0, then dH]) and ([49]) become 



V(26 + cr s ) + 3cr i 



B' 
AB 



0. 



Vo s + a s ^ = 0. 



Integrating ([55]) we obtain 



B 



(82) 



(83) 



(84) 



where f(t) is an arbitrary integration function. However 
since we have cr s (t, 0) = we must put f(t) = implying 
that a s = too. Also, from (gg) VP = 0. 

Inversely, if the fluid is shearfree and geodesic then it is 
necessarily nondissipative. Indeed, if the fluid is shearfree 
this condition can be integrated (in general, not only in 
the geodesic case) to obtain (see [14j for details) 



B = Ab(r) and C = Ac(r). 
where b and c arc arbitrary functions of r. 



(85) 
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Then, feeding back the above expressions for B and 
C into (| 13[) and taking into account ([76")) it follows that 
the expansion scalar is homogeneous (V0 = 0), implying 
because of ([4*8")) that geodesic shearfree fluid is necessarily 
non-dissipative . 

Thus we have proved that for the geodesic fluid, q a = 
H a p = (T Q( g = 0. Whereas in the general case we 



have stablishcd that a 







V6 — nq = H, 



a/3 



0. 



A' / B' 




B' 


a 


A \B 


+ c) 


+ B 


c 


A" 


C" 




2 


-A + 


~c ~ 







= kP z A\ 



(87) 



VIII. THE INHOMOGENEITY FACTOR AND 
ITS EVOLUTION 

In the spherically symmetric case it has been shown 
that in the absence of dissipation the necessary and suf- 
ficient condition for the vanishing of the (invariantly de- 
fined) spatial derivative of the energy density is the van- 
ishing of the scalar associated to the trace free part of 
X a p . For obvious reasons such a quantity was called the 
inhomogeneity factor. In other words, the inhomogeneity 
factor (say is that combination of physical and geo- 
metric variables, such that its vanishing is a necessary 
and sufficient condition for the homogeneity of energy 
density (if dissipation is present then additional terms 
including dissipative flux appear (57l|V 

The extension of such a definition to situations where 
there is only one relevant spatial coordinate (as in the 
case considered here) is rather obvious: 

f = V/i = 0. 

The two equations ([50]) and ([51"]) determine the inho- 
mogeneity factor, whereas ((52)) and ([53"]) describes its evo- 
lution. 

Nevertheless, in the present case the situation is much 
more complicated (than in the spherically symmetric 
case) due to the fact that not only the tensor X a /3 is 
expressed in terms of three structure scalars (instead of 
two in the spherically symmetric case) but also due to the 
fact that the magnetic part of Weyl tensor is in general 
not vanishing and appears as a possible source of energy 
density inhomogeneity. 

Because of the above reasons, we were unable (in the 
general case) to identify explicitly the inhomogeneity fac- 
tor. The very simplified cases where this was possible are 
not very interesting and therefore we shall not include 
them here. 



IX. ALL STATIC ANISOTROPIC CYLINDERS 

We shall show in this section that all possible solutions 
of the static case are completely determined by a triplet 
of structure scalars. This is a reminiscence of the spher- 
ically symmetric case, where all possible static solutions 
are determined by a couple of structure scalars [l[ . 

In the static case the field equations read: 



A" B" 



A' 



A • B U 1 = K ^ A 



Then introducing the auxiliary variables 
A' B' C 

e ~s' C "C' 

we can write the field equations as: 

-e - e - c - c 2 + + - cc = ^a 2 , 



uj£ + + £C = nP r A 2 



or. 



u' + C + C = *PzA*, 



(89) 



(90) 



(91) 



(92) 



(93) 



(94) 



(95) 



+ e -oj(;-oj(~tt = K(P^-Pr)A 2 , (96) 



to' + C + C - w£ - u>£ - £C = k{P z - P r )A\ (97) 



e+e-c -c 2 = k(p^-p z )a 2 . 

The scalars E s y Ek take the form 



1 

2A 2 
I 

2A 2 



E s = ^{-Lo' + C + C 2 +^-uj(-tiC), (99) 



E k = — (-^' + £' + £ 2 -^ + ^C-£C) ■ (ioo) 



Then, using ([32]), ([33]), ([M]), (JET]), flU and ((M)]) we 

can write 



B" 


C" 


A 1 




C'\ 


B' 


a 


~~B 


~ ~C ' 


h A 


(§- 


h c) 


B 


c 



= KfiA 2 , (86) 



Y S A 2 = - J + w£, 
Y k A 2 = - J + u>(. 



(101) 
(102) 
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Integrating (|101[) wc obtain 
A = a cxp / B 



-Y S A 2 
B 



where a is a constant. Thus for any given Y s , we obtain 
from (|103p a relationship between A and B, 



dr. 


(103) 


A 1 ( B' 




b' a 






A \ B 


' c) 


+ B C 



B = B(A) ^>uj = 



Next, from ([TU21 

A = 7 cxp / C 



-Y k A 2 
C 



dr I dr, 



(104) 



(105) 



where 7 is a constant. Therefore, giving Y k we obtain a 
relationship between A and C, 



C = C(A) =>u = u(Q. 



(106) 



Then from (|104|) and (|106|) we can express any of (ui, £, () 
in terms of the other two. Therefore in (|98[) we can ex- 
press £ and £ in terms of to. obtaining a differential equa- 
tion for uj. This can be solved for a given P^ — P z and 
once ui is obtained, we can get £ and £ from (| 104[) and 

(Hog). 

Alternatively we may use (|96[) or (|97|) in which case we 
should provide either — P r or P z — P r . Once the metric 
functions are found we can obtain physical variables from 
field equations. 

Now, from (f33l) and (EZl) it follows 



(107) 
(108) 



K(Pt-P r ) = -(Y k + X k ), 
and from ([35]) and (|5H]| 

K(P z -P r ) = -(Y s + X s ). 



Therefore any static anisotropic solution is determined 
by a triplet of scalars (Y k ,Y s ,X k ) or (Y k ,Y s ,X s ). 
We shall next consider some special cases. 



AT_ 
~A 

AT_ 
'IT 



~~C 

El 
~B 



.4 



= kPA\ 
= kPA 2 , 



7T ' 



(109) 
(110) 
(111) 



Then, from (|110[) and (|lllj) we obtain 
C" B" 



(112) 



which in terms of the auxialiary variables introduced be- 
fore, read 



which is a Ricatti equation for £ ( or £). 

The general solution of (|113|) takes the form 



(113) 



1 



k{r)' 



(114) 



with 



k(r) = e 2 ^ dr f^j e~ 2 ^ dr dr + aj . (115) 



Turning back to B and C we can integrate (|114[) , to ob- 
tain 



C = PBe 



tlr 



(116) 



where f3 is a constant. From regularity conditions we 
must impose that C(t, 0) = for any acceptable solution. 



A. Isotropic cylinders 

In this case P r = P z = P^ = P and field equations 
(E7l)-([89l) become: 



B. Conformally flat solutions 

In the static case we obtain from (l24l) 



El 
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\ A 'l 


f B< 


A 1 


C'\ 


B' C 


A" 


C"~ 




A \ 
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B C 


-a' 


~c ' _ 
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'A' 




A 1 


C'\ 


b' a 


A" 


B" 


~2A2 


A 
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+ B C 


f A 


B 



(117) 



(118) 



Then, defining two functions c(r) and b(r) such that C = 
Ac(r) and B = Ab(r) the conformal flatness condition 



(E s = E k = 0) implies 



(119) 
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and 

b" = -b'. (120) 

c 

From (|120|) wc get d — r yb' (where 7 is a constant) and 
feeding this back into (|119j) we obtain 

b(r) = /3cos(er), (121) 

where e and /3 are constants, and the regularity condition 
c(0) = has been used. 

In the isotropic case we obtain from ([1091 lllll) 

A = constant; n = — 3P = constant. (122) 

All these, conformally flat, solutions have been de- 
scribed in detail in [58j , including the discussion on junc- 
tion (Darmois) conditions on the boundary surface. 

X. DISCUSSION AND SUMMARY OF RESULTS 

A comprehensive study on cylindrically symmetric rel- 
ativistic fluids by means of structure scalars have been 
carried out. 

We have first defined the complete set of such scalars 
corresponding to our problem. It turns out that there are 
eight structure scalars (Xt, Xk, X s , Yr, Y s , Y/~, Z q , Zh) in 
contrast with the spherically symmetric case where there 
are only five. Besides, two scalars defining the shear ten- 
sor ((Tfc,(T s ) and three scalars defining the electric and 
magnetic parts of Weyl tensor (E)~, E s , H) were also in- 
troduced. 

Next we have stablished a set of equations governing 
the structure and evolution of the system under consid- 
eration and brought out the role of structure scalars in 
those equations, in order to exhibit the physical relevance 
of the former. 

We have first considered the dynamical equation (fP| 
derived from conservation laws. We have next coupled 
the above mentioned equation with a transport equation 
derived from a causal dissipative therory. The resulting 
equation exhibits the decreasing of the effective inertial 
mass term due to thermal effects. 

Next, we have brought out the relevance of structure 
scalars. Our main results in this respect can be summa- 
rized as follows: 

• Three of them Z q , Zh and Xt, have an evident 
physical meaning, and therefore do not require fur- 
ther discussion. 

• Yt has been shown to control the evolution of the 
expansion scalar through the Raychaudhuri equa- 
tion P5]) . whereas Y s and Yfc control the evolution 
of the shear through (gBJ) and (|47|) . 

• A very tight link between the shearfrec condition, 
Z q and Zh appears from equations (|15j) and (flUl) . 
Thus, it has been shown in the geodesic case that 



necessary and sufficient conditions for the fluid to 
be shearfree are Z q = Zh = 0. In the general case 
it has been shown that the shearfree condition im- 
plies the vanishing of the magnetic part of the Weyl 
tensor and a direct relationship between the inho- 
mogencity of the expansion scalar and the dissipa- 
tive flux. The former result explains the absence of 
gravitational radiation in the shearfree case already 
commented in (T3 |. 

• The two equations (|50"]l and (|5Tj) relates X s , A% Y B , 
Yfc with energy density inhomogeneity and there- 
fore should provide a definition of the inhomogene- 
ity factor(s) in terms of some structure scalars. In 
the same order of ideas equations ((52)) and (|55|) 
describe the evolution of such factor (s). Unfortu- 
nately, in the general case, we were unable to iso- 
late them (it was possible only in very simplified 
situations) . 

• We have next considered the static case. The 
main result from this section is the obtained proce- 
dure allowing to determine any possible solution in 
terms of a triplet of structure scalars (Yk, Y s ,Xk) or 
(Yk,Y s ,X s ). Particular subcases such as isotropic 
or conformally flat cylinders were studied, obtain- 
ing specific restrictions about the existence of solu- 
tions satisfying regular conditions on the symmetry 
axis. 

As expected, in a general study as the one presented 
here, a great deal of questions remains unanswered. Thus 
before ending we would like to present a partial list of 
issues that should be addressed in the future: 

• From (fi2"j) it follows that the "gravitational" term 
vanishes not only if H = but also if Ek = E s . 
Why? What else does this latter condition imply? 

• We have seen that the shearfree condition implies 
the vanishing of the magnetic part of the Weyl ten- 
sor. What are the implications on the shear of the 
vanishing of the magnetic part of the Weyl tensor 
in the general case? 

• Are there purely magnetic solutions? 

• Could it be possible to find the exact solution cor- 
responding to nondissipative dust with shear (the 
analog of the Lemaitre-Tolman-Bondi solution)? 
Would this solution have a nonvanishing magnetic 
part of Weyl tensor? 

• We have identified the subset of equations which 
should determine the inhomogeneity factor and its 
evolution, but we were unable to isolate such a fac- 
tor in the general case. Is this possible? 
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APPENDIX I 



The nonzero components of the Einstein equations 
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APPENDIX II 



G a p = nT a p, 



The non null components of the Weyl tensor Cap-ys for 
(123) © are 
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The two components of the Riemann tensor appearing in 
(El) are 
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2323 



/ .BC 



BC B'C 
BC ~ ~B~C 



(134) 
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XI. APPENDIX III 

From their definition it is not difficult to express struc- 
ture scalars through metric functions and their deriva- 
tives, these expressionas are: 



1 AB AC BC B" C" 
Xt ~ A* \AB + AC + BC ~ ~B ~ ~~C 



A'B' A'C B'C 
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AC BC 
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